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In thermal phases, the quantum coherence of individual degrees of freedom is rapidly lost to the 
environment. Many-body localized (MBL) phases limit the spread of this coherence and appear 
promising for quantum information applications. However, such applications require not just long 
coherence times but also a means to transport and manipulate information. We demonstrate that 
this can be done in a one dimensional model of interacting spins at infinite temperature. Our 
protocol utilizes protected qubits which emerge at the boundary between topological and trivial 
phases. State transfer occurs via dynamic shifts of this boundary and is shown to preserve quantum 
information. As an example, we discuss the implementation of a universal, two-qubit gate based 
upon MBL-protected quantum state transfer. 
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The faithful storage and transmission of quantum in¬ 
formation is the basic building block for applications 
ranging from information processing to communication 
and metrology [1-3]. In a traditional architecture, quan¬ 
tum bits connect via channels that coherently shuttle 
information between remote nodes [4]. Constructing 
such a channel from an interacting many-body system 
at high temperature is generally considered impossible; 
once quantum information disperses into the system, it 
rapidly decoheres due to scattering with thermal excita¬ 
tions. To avoid this, typical quantum channels, such as 
mechanical resonators [5, 6], optical photons [7, 8], su¬ 
perconducting strip-lines [9, 10], and spin chains [11-26], 
are either specially tailored few-body systems or operate 
at ultra-low temperatures in order to freeze out parasitic 
degrees of freedom. 

Recently, an alternate strategy to deal with thermal 
excitations has emerged: localization [27-31]. The intro¬ 
duction of sufficiently strong quenched disorder leads to 
a many-body localized (MBL) phase [27-62] with effective 
Hamiltonian, 

H = J2Jir,^ + Yl ^ + ■ ■ • ( 1 ) 

i ij ijk 

where rf are local integrals of motion and ffij... de¬ 
cay exponentially [32-35]. In general, these rf are lo¬ 
cally dressed versions of the underlying physical degrees 
of freedom. From Eqn. (1), it is clear that their 2 :- 
polarization never decays (i.e. Ti ^ oc). However, su¬ 
perpositions dephase rapidly due to the Ising-type inter¬ 
actions with neighbors (i.e. T 2 ^ 1/*^)- principle, 
for each rf this dephasing can be canceled by a local 
spin-echo (i.e. T 2 ^ oc) [35, 39]. In practice, there are 
three significant challenges: 1) there is only a finite ran¬ 
dom overlap between the physically addressable degrees 
of freedom and the rf, 2) quantum information is only 
protected in the presence of continuous spin-echo, 3) it 
is unclear how to achieve robust quantum state transfer 
and logic gates. 



FIG. 1. (a) Schematic representation of the boundary be¬ 

tween a trivial and SPT insulator. A single qubit is bound at 
the interface. Strong disorder (black line) leads to many-body 
localization in both phases, (b) Shuttling of the edge qubit 
occurs via shifts in the trivial-SPT boundary. 


In this Letter, we demonstrate the storage and trans¬ 
mission of quantum information in an interacting, many- 
body spin chain at infinite temperature (i.e. finite energy 
density) [64]. Our approach requires neither initialization 
nor continuous spin-echo. Rather, we exploit the bound¬ 
ary qubit at the interface between a many-body localized 
trivial and symmetry protected topological (SPT) phase 
(Fig. 1) [35-38]. Symmetry prevents the operators corre¬ 
sponding to this degree of freedom from appearing in the 
effective Hamiltonian [Eqn. (1)], implying that T| ^ oc 
even in the absence of spin-echo. Our main result is sum¬ 
marized in Table I (bottom row); the 0(1) local spec¬ 
tral gap of the MBL phase is sufficient to protect static 
storage (Eig. 2a), dynamic manipulations (Eig. 2b), and 
quantum gates (Eig. 2c) between boundary qubits; we 
emphasize that our approach does not rely upon, nor 
work in, the exponentially slow, strictly adiabatic limit 
[63]. 

Protocols and model —In one dimension, the canonical 
example of an SPT phase is provided by the Haldane 
phase, whose boundary binds a single qubit (Eig. la) 
[66, 67]. But what exactly does this mean? In one sense, 
this means that coherent quantum information may be 
stored and retrieved from the edge qubit. Let us de¬ 
fine a “static” protocol to sharpen this: measure the a- 
component, of the qubit to initialize it, wait for a 
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TABLE L MBL protected storage (static) and shuttling (dy¬ 
namic) 



Static 

Dynamic 

Zero T, clean 

/ 

/ 

Finite T, clean, non-interacting 

/ 

X 

Finite T, clean, interacting 

X 

X 

Finite T, disordered, interacting 

/ 

/ 


time r, and then measure again. At zero tempera¬ 
ture, the measurements coincide and there is neither loss 
of polarization nor coherence. 

Being bound to the edge also implies that the qubit 
moves when the boundary shifts (Fig. lb). But, does this 
motion decohere the quantum information? Let us define 
a corresponding “dynamic” protocol: measure move 
the boundary of the SPT phase (e.g. by ramping local 
fields and converting a region into the opposite phase) 
on a time-scale tr^ hold for a time r, ramp back and 
measure. At zero temperature, the measurements again 
coincide when tr ^ 1/A^, where A is the bulk gap, as 
required by the quantum adiabatic theorem (Table I, top 
row). 

The introduction of finite temperature invalidates the 
above lore. Surprisingly, a distinction emerges between 
the ‘static’ and ‘dynamic’ protocols (Table I). Consider 
the following microscopic spin Hamiltonian [36], 

h = Y, \i<yU<yi<yUi + E + E ( 2 ) 

i i i 

where cr" are Pauli operators. This model has a Z 2 x Z 2 
symmetry corresponding to products of on the even 
and odd sites. The symmetry protects the Haldane phase 
realized when A is the dominant coupling. A trivial in¬ 
sulator arises when h ^ X,V and thus a trivial-SPT 
boundary can be realized (Fig. la) by spatially varying 
the ratio X/h. In the simplest case when V = 0 one 
can write down an explicit expression for the edge qubit. 
Setting hi = 0 on the SPT side (i > 1) and = 0 on 
the trivial side (i < —1), the Hamiltonian is a sum of 
commuting terms. One finds that the following Pauli op¬ 
erators at the boundary commute with the Hamiltonian, 
E^ = (Jocrf, E^ = (Tgcrf, E^ = cTo? t)ut are odd under the 
Z 2 X Z 2 symmetry and hence do not appear in H. The 
Pauli algebra guarantees that they generate a two fold de¬ 
generate edge state. For a more generic Hamiltonian, the 
edge states are dressed by exponentially localized tails. 

To analyze the protocols, we switch to a convenient 
representation of the Hamiltonian in terms of Majorana 
fermions: 

= E iXaf-iil+i + E Valihi+iXi+i- 

i i i 

( 3 ) 
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FIG. 2. Exact diagonalization and numerical integration of 
Hamiltonian evolution in the interacting, disordered, infinite 
temperature case. All numerics are performed using a box 
distribution with mean and width shown in the figure, (a) 
Depicts the static coherence time of the edge qubit, which 
grows exponentially with system size. The total system size 
is L = ^ + 2 with a trivial insulator (A = 0) occupying the 
first two sites followed by an SPT region of length i. (b) 
Fidelity of the dynamic protocol. The boundary qubit is suc¬ 
cessively shuttled inwards by a total of four lattice sites with 
the ramp profile shown in the bottom panel, (c) Shows 145 
disorder realizations of the dynamic protocol for L = 16. In 
each realization, the boundary qubit is shuttled by a total of 
eight lattice sites in four successive steps. The pSWAP gate 
time improves exponentially while the adiabatic ramp time is 
governed by the 0(1) local spectral gap. 


7 ^ and 7 ^, with the Pauli algebra given by, af = 

= Y{j<i (G^i) = Y{j<i (iij'yj) il- 

Finite temperature, elean —We begin the finite tem¬ 
perature analysis with H = 0 in the clean case (uniform 
couplings), where the Hamiltonian is quadratic and has 
a simple pictorial representation (Fig. 3a). This clari¬ 
fies the nature of the Z 2 x Z 2 symmetry: the fermions 
hop on two chains, each of whose parity is conserved. 
In the fermionic language, the edge qubit operators are 
given by the zero modes ( 70 , 7 i triangles. Fig. 3a) along 
with Jordan-Wigner strings that lie entirely in the trivial 
phase [65], 


s- = = n (nhl) 

j<0 

= a^at = -no7i 

^^ = ^^0 = 11 Xo- (4) 

j<0 


Here, each spin is represented by a pair of real fermions 


These strings reflect the SPT nature of the qubit; they 
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attach to all local operators odd under the Z 2 x Z 2 sym¬ 
metry. 

In the static protocol, both the zero modes and the 
Wigner string are conserved and thus coherence persists 
even at finite temperature (Table I, second row). The 
dynamic protocol is more subtle. To shift the bound¬ 
ary, one ramps up Ao(t),A_i(t) as shown in Fig. 3c. At 
all times during this ramp, there is a unique zero mode 
with an 0(1) gap to the bulk for each chain; thus, if each 
ramp proceeds slower than this gap scale, the zero mode 
simply shifts and returns with no phase accumulation. 
However, the last four fermionic operators of the Wigner 
string ( 7 ^ 257 - 257 - 157 - 1 ? 3a) are carried across the 

phase boundary and necessarily move into modes above 
the gap. At T = 0, there are no fermionic excitations in 
the system that can be carried across with these modes; 
however, at finite temperature, any excitation carried 
across the boundary disperses into the bulk of the SPT 
phase. Thus, the qubit’s coherence is preserved through 
the dynamic protocol only at T = 0 in the clean system. 

To confirm this picture, we consider the trace fidelity 
of the edge qubit, 

F°‘{t) = ^Tr [e“'®"S“(t)S“(0)] , (5) 

Zj 

as a function of time, where Z is the partition function 
and p = l/ksT. This fidelity is precisely the expected 
correlation of the measurements in the static and dy¬ 
namic protocols. We numerically simulate the dynamic 
protocol at infinite temperature on finite size systems of 
length L with the boundary at L/2. For large systems, 
we expect the dynamic protocol to fail (Fig. 3b), while for 
sufficiently small systems, it is possible for the ramp rate 
to resolve the mini-gaps ^ \jL within the bulk bands. 
Indeed, for fixed the qubit’s coherence is lost as L 
increases (Fig. 3b inset) as neither YZ nor recover 
(Table I, second row). The persistence of reflects its 
lack of a Wigner string and further confirms our above 
analysis. 

In the interacting case, V ^ 0, the situation is much 
worse: both the static and dynamic protocols fail at finite 
temperature (Table I, third row). Indeed, scattering by 
bulk excitations leads to a static coherence time T 2 ^ 
\/V^ [71]. Until recently, this would be the end of the 
story. However, it is now believed that strong disorder, 
leading to many-body localization can protect SPT order 
even at finite temperature [35-38]. 

Finite temperature, disordered —The localization of 
bulk excitations suppresses scattering and enables the 
static coherence of the edge qubit to persist in this MBL 
SPT phase. More specifically, in the trivial insulator, 
the local integrals of motion are given by rf = UafU^, 
while in the SPT insulator, rf = where 

1/ is a local Z 2 x Z 2 symmetric unitary that diagonalizes 
the Hamiltonian. These symmetry-even local operators 
appear directly in Eqn. (1). As previously discussed, un- 
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FIG. 3. (a) The U = 0 limit in which H reduces to a pair of 
decoupled fermion chains. The SPT insulator (on the right, 
z > 0) has A (black line) h (dotted line) while the trivial in¬ 
sulator (on the left, z < 0) has A (gray line) h (dotted line). 
Ramping the boundary terms Ao,-i, shifts the interface and 
hence the location of the edge qubit, (b) Trace fidelity of the 
dynamic protocol in the infinite temperature, non-interacting, 
clean case for a chain of length L = 40 with the boundary at 
L/2. On the SPT side, /z = 0.1, A = 1 while on the trivial 
side, h = 0.1, A = 0.01. The couplings Ao and A_i ramp 
linearly from 0.01 to 1 as a function of time [dotted line in 
(c)]. (inset) For L = 40, a ramp time of tr ~ 10^ is slow 
enough to resolve the mini-gaps ~ 1/L within the band; thus, 
the dynamic protocol exhibits a nearly perfect fidelity. How¬ 
ever, this is a finite size effect and increasing the system size 
immediately leads to a failure of the protocol, (c) Same as 
(b) but with strong disorder, ^A = 0.7A, Sh = h. Localization 
of the Wigner strings enables a recovery fidelity > 99% with 
the same ramp time. 


like the bulk rf, the boundary qubit, = UY^W is 
symmetry odd and therefore cannot appear in Eqn. (1). 
Instead, the boundary qubit can only be dephased by 
interactions involving other edge modes. Eor an SPT 
region of length I, this implies T| ^ e^, as confirmed 
numerically in Eig. 2a. 

But, can localization protect coherence during dynam¬ 
ical protocols? The non-interacting case already provides 
a tantalizing hint that this may be. Eigure 3c shows the 
successful recovery of quantum information in a disor¬ 
dered system at infinite temperature. We note that the 
ramp time used is identical to that of the failed attempt 
in the clean system (Eig. 3b). Here, localization plays 
a crucial, but different role than in the static case. It 
prevents the dispersion of the piece of the Wigner string 
carried across the boundary, limiting its spread to a lo¬ 
calization length ^ independent of system size. Thus, so 
long as the ramp time is adiabatic with respect to the 
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local gap, tr ^ 1/A^, the string returns to its original 
position following the protocol. 

In the interacting case, one expects that the adiabatic- 
ity condition is modified by at most a factor 2^ ^ 0(1). 
Numerics confirm this expectation: The top panel of 
Fig. 2b depicts the fidelity of a successful dynamic proto¬ 
col at infinite temperature, with both strong interactions 
and disorder. Two successive ramps (ramp profile shown 
in Fig. 2b bottom) move the edge qubit a total of four lat¬ 
tice sites before holding and returning; the system size is 
L = 12 with the trivial phase initially occupying the first 
two sites. The example shown is typical among ^10^ 
samples; failure (less than 90% recover fidelity) arises for 
less than 3% of the total disorder realizations and results 
mainly from two causes: 1) when there are additional, 
unexpected, trivial-SPT boundaries and 2) when the lo¬ 
calization length is significantly larger than a few lattice 
sites. 

MBL protected quantum gates —Between the static and 
dynamic protocols, we now have robust control over a 
single edge qubit in the MBL SPT phase. This raises 
the question of whether one can perform a non-trivial 
gate between two such qubits. For a finite SPT region 
of length the two edge qubits interact with each other 
through their exponentially localized tails. In the clean, 
non-interacting case, this leads to coherent oscillations 
on a time scale, tq = f (^)^^^~^ [71]. 

The nature of these oscillations is apparent from the 
fermionic Hamiltonian in Eqn. (3) at H = 0. At time 
To, the zero modes at each end of the SPT region have 
swapped position. In the language of the edge qubits, 

^l(^) ~ cos(0/^t)S|^(O) -|- sin(0/^t)iPS^(O) 

T.l(t) = cos((/>Bt)SnO) - sin((/>Bt)iPS|,(0) 

= - sin(0Bt)iPS|(O) + cos(0Bt)E|^(O) 

T.R(t) = sin(0A^)^PSi(O) + cos(0At)S^(O), (6) 

where P,P index the left and right qubits, 0 a, 0b are 
the oscillation frequencies (0 = 27r/To for the clean case), 
and P = ]\i^i parity operator for the SPT region. 

At time tq, Eqn. (6) produces a non-trivial gate be¬ 
tween the edge qubits, which we term: pSWAP. It cor¬ 
responds to a SWAP gate between the edge qubits en¬ 
tangled with the parity of the SPT region. This entan¬ 
glement is rather benign since P^ = I, implying that a 
second pSWAP cancels the unwanted entanglement. Un¬ 
fortunately, the naive pSWAP gate suffers from two is¬ 
sues: 1) it does not directly lead to universal quantum 
computation and 2) it is exponentially slow in the sepa¬ 
ration of the edge qubits. 

The first issue is easily remedied by placing an ancilla 
qubit next to each SPT edge qubit. With full control over 
this quantum register, we can perform an MBL-protected 
universal quantum gate, as shown in the quantum circuit 
of Eig. 4. The explicit Heisenberg evolution of the oper¬ 
ator basis is indicated below the circuit (P = 6). Erom 
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FIG. 4. Quantum circuit demonstrating a universal MBL- 
protected gate between ancillae qubits LM and RM. The 
Heisenberg evolution of all operators of an L = 6 SPT region 
is shown following each gate. At the end of the evolution, 
a universal controlled-T gate has been implemented between 
LM (control) and RM (target). 


this, one can deduce that a universal controlled-Y gate 
occurs between the two ancillae. 

The second issue can be remedied by using the dy¬ 
namic protocol for quantum state transfer. By shuttling 
the edge qubits closer to one another, the pSWAP gate 
time is exponentially improved. The coherence time of 
the SPT edge qubits ^ e^. Meanwhile, each shuttling 
step takes time ^ 1/h^; following m steps, the pSWAP 
gate occurs in a time ^ giving an exponential en¬ 

hancement between gate and coherence times. This is 
explicitly shown in Fig. 2c; we simulate 145 disorder re¬ 
alizations on an interacting L = 16 chain at infinite tem¬ 
perature. We shuttle the left SPT qubit from 2 to 8 
sites. At each separation, we extract the pSWAP time 
and scatter this against the total ramp time. Each ramp 
is optimized to ensure a final fidelity F > 90%. In the dis¬ 
ordered case, since 0 a % 0b in general, one must match 
the phase accumulated by and [Eqn. (6)]. This 
is accomplished by adjusting the hold time between each 
ramp. 

Conclusion —The binding of a qubit to the interface 
between trivial and SPT insulators has enraptured a gen¬ 
eration of theoretical physicists. At zero temperature, 
shifting the position of this interface coherently drags 
the qubit along, analogous to the shuttling of Majorana 
fermion zero modes in topological superconducting wires 
[68, 69]. In clean systems at finite temperature, neither 
storage nor shuttling of the SPT qubit works. We have 
demonstrated that strong disorder, leading to many-body 
localization, protects qubit coherence in both “static” 
and “dynamic” protocols. Localization plays two roles: it 
prevents scattering from thermal excitations in the bulk 
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and also prevents the “Cherenkov” emission of quantum 
information during transport. This enables the coher¬ 
ent control of qubits in disordered, strongly interacting 
systems at high temperature. Furthermore, our con¬ 
struction of the pSWAP gate demonstrates that these 
qubits may be wired up and universal gates performed 
between them. Looking forward, it would be interesting 
to consider a wire network of disordered SPT insulators 
as a platform architecture [68]. The most natural system 
where our approach offers a competitive advantage is one 
in which Hamiltonian engineering is possible but cooling 
into fiducial initial states is difficult. One intriguing pos¬ 
sibility is in driven systems where “stroboscopic” SPT 
phases can emerge upon periodic modulation of interac¬ 
tions [70]. In these cases, many-body localization is, by 
itself, unable to provide coherent local degrees of free¬ 
dom, thus motivating the need for an MBL SPT qubit. 
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Supplemental Material for Many-body Localization Protected Quantum State Transfer 


N. Y. Yao, C. R. Laumann, A. Vishwanath 


Static coherence times—Here, we provide additional numerics and analytic details regarding the state transfer protocol. First, 
we consider the coherence time of the edge qubit in the infinite temperature, clean, interacting system. Fig. Sla depicts the 
trace fidelity as a function of time for a length L = S SPT chain. Fitting the decay envelope, , of allows one to 

extract T^. Moreover, by varying the interaction strength V, while holding all other parameters fixed, we extract the scaling, 
T 2 ^ X/V‘^ (Fig. Sib). Adding strong disorder brings us to the case of the maintext, where T 2 ^ e^; an example of this 
extended coherence is shown in the inset of Fig. Sla. 




FIG. SI: Infinite temperature, clean, interacting case, (a) Quantum information storage in the SPT edge qubit as a function of time with L = 8, 
A = 1, /i = 0.05, V = 0.01. (inset) Same but for uniform disorder of width SX = 0.7, Sh = 0.05, SV = 0.01. (b) Inverse quadratic scaling 
of the static coherence time as a function of interaction strength V. 

Non-adiabatic transitions—To further confirm that the dynamic protocol fails in the clean case owing to dispersion of the 
Wigner string into the bulk, we consider an isolated (non-interacting) SPT chain, as shown in the inset of Fig. S2. Without the 
semi-infinite trivial insulator to the left, the edge qubit Pauli operators are simply given by, 

= o-Q = 7o 
El = 

Sf =^i. (SI) 

To perform an analogous dynamic protocol, we shuttle the edge qubit inwards by ramping down the first two topological terms 
Ao (Fig. S2). We hold for a time rw = 100/A and then ramp Aq back to its original value. As in the maintext, at all times during 
the protocol, there is a unique zero mode with an 0(1) gap (^ h) to the bulk for each chain; thus, if each ramp proceeds slower 
than this gap scale, the zero mode simply shifts and returns with no phase accumulation. Moreover, since the edge operators no 
longer carry Wigner strings, this is the only relevant adiabaticity condition. 

The trace fidelity of this static protocol is shown as a function of the ramp time in Figure S2. The system size is L = 20 and 
A = 1, h = 0.05. One finds that the functional dependence on the ramp time is well characterized by a Landau-Zener-like 
formula. 


F w 1 - (S2) 

with 77 fs 2 . 877 . Note that the fidelity of this “no Wigner string” dynamic protocol is greater than 95% for all tr > 150/A. This 
is in stark contrast to Fig. 3a of the maintext where all parameters are identical except for the existence of a finite trivial chain to 
the left of the SPT region. The existence of this trivial chain implies the the edge qubit operators carry Wigner strings and thus 
require a significantly more stringent adiabaticity condition that resolves the bulk mini-gaps. 

MBL-protected universal two-qubit gate—We derive the pSWAP enabled universal two-qubit gate between the SPT edge 
ancillae. As in the maintext, let us assume full control of an edge quantum register that contains the MBL SPT edge qubit and 
an ancillae qubit. The quantum circuit shown in Fig. 4 of the maintext is calculated for the non-interacting (V = 0) case by 
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FIG. S2: Fidelity of the dynamic shuttling protocol for a non-interacting SPT chain as a function of the ramp time, (inset) Pictorial represen¬ 
tation of the Majorana chains; note that there is no boundary to a trivial insulator. Also shown is the ramp profile for Aq. 


following the Heisenberg evolution of the modes. For simplicity, we neglect the dressing of the edge operators by exponential 
tails. At the end of the circuit, the ancillae operators have transformed as: Xl XlZr, Zl ^ ZlZr, Xr —YlXr, 
Zr Zr. Crucially, they are now disentangled from the intermediate SPT chain. From the evolution, we can analytically 
compute the underlying unitary gate [1]. Ultimately, this guides the intuition, but the precise gate that is implemented in the 
interacting case cannot be derived in this fashion; rather, our full numerics demonstrate a universal entangling gate with a finite 
fidelity reduction for U 7 ^ 0 . 


For the non-interacting case, we explicitly derive the action of the unitary on the basis {tt? it? tt? it} (of loft and right ancillae 
qubits). Let us start with |tt)- This is a +1 eigenstate of both Zr and Zr and thus, it must also be a +1 eigenstate of the final 
transformed operators ZrZr and Zr. This immediately implies that U |tt) = ^ Itt)- The other three cases are given by. 


|;t) = in) 
u in) = uxlU^u in) 
u in) = xlZhu in) 
u in) = aXiZR in) 
u in) = «in) 

in) = in) 
u in) = uxru^u in) 
u in) = -ylXru in) 
U in) = -ia in) 

in) = in) 

U\ii) = UXLU^U\n) 
u\ii) = XLZRU\n) 

U\ii) = -iaXLZR\i\) 

U\U.)=ia\n)- 

Thus, the final unitary gate between the ancillae qubits is given by, 

/I 0 0 0 \ 

I 0 1 0 0 I 
^ ~ 0 0 0 • 
Vo 0 i 0/ 


(S3) 


(S4) 


(S5) 


(S 6 ) 
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which is a controlled-F gate where the control is right ancillae and the target is the left ancillae qubit. This gate is universal. 


[1] D. Gottesman, Proceedings of the XXII International Colloquium on Group Theoretical Methods in Physics, eds. S. P. Corney, R. Del- 
bourgo, and P. D. Jarvis, pp. 32-43 (1999). 



